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Abstract 

Peptide Nucleic Acids (PNA) are non-natural oligonucleotides mimics, 
wherein the phosphoribose backbone has been replaced by a peptidic moi¬ 
ety (N-(2-aminoethyl)glycine). This peptidic backbone lends itself to sub- 
stitntion and the y-position has proven to yield oligomers with enhanced 
hybridization properties. In this study, we use Nuclear Magnetic Reso¬ 
nance (NMR) and Circular Dichroism (CD) to explore the properties of the 
supramolecular duplexes formed by these species. We show that standard 
Watson-Crick base pair as well as non-standard ones are formed in solution. 
The duplexes thus formed present marked melting transition temperatures 
substantially higher than their nncleic acid homologs. Moreover, the presence 
of a chiral gronp on the y-peptidic backbone increases further this transition 
temperature, leading to very stable duplexes. 

PNA duplexes with a chiral backbone present a marked chiral secondary 
structure, observed by CD, and showing a common folding pattern for all 
studied structures. Nevertheless small differences are observed depending on 
the details of the nncleobase sequence. 
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Introduction 


Peptide Nucleic Acid (PNA) are synthetic oligonucleotides first reported 
by Nielsen[T] with a backbone that recapitulate DNA’s inter-nucleobase dis¬ 
tances. PNAs were developed as DNA mimics in order to recognize DNA 
double helix [ 2 ] and to form hybrid duplexes with DNA[3j or homoduplexes jl]. 
PNAs duplexes have high thermal stability, are not sensitive to nucleases or 
protases, and are metabolically stable. However, as polyamide backbone has 
no charge, no Coulomb repulsion is reported and PNA aggregates may exists 
in solution, and their solubility in water is limited|3]. PNAs have attracted 
interest in molecular biology and numerous applications have been reported, 
we can cite beacons for duplex DNA[ 6 ], and PNA as a diagnostic tool[7j, a 
genomic tool[H], or a supramolecular barcodes [UlHD]- 

These applications were leveraged on the unique properties of PNA hy¬ 
bridizing to DNA or RNA. The enhanced hybridization properties of 7 - 
modified PNA have been recently harnessed for gene editing apphcations|llj. 
PNA have also been used to tag small or macromolecules and program their 
assemblies based on hybridization j9] [TO]. PNA homoduplexes have also been 
used in programmed assemblies as recently illustrated for the programmed 
pairing of PNA-tagged protein fragments [T^. While interactions of PNA 
with DNA or RNA have been extensively studied, homoduplex formation of 
PNA are not well characterized, particularly for modified PNAs. 

PNA-DNA hybridization by Watson-Crick base pairing are well studied[3] . 
They reveal a helix formed by the 7 -carbon and the nitrogen of the ter¬ 
tiary amide of the PNA backbone [TH]. stabilized by the sequential base 
stackinglll]. PNA neutral backbone bring thermal stability in PNA/DNA 
duplex compared to DNA/DNA duplexes, and present a better specificity [I5]. 
Because of the absence of charge on the PNA backbone, the PNA/DNA hy¬ 
brids are not dependent on the ionic strength of the solvent in contrast to 
DNA homoduplexes. Because of a flexible backbone, Watson-Crick, Hoog- 
steen, reverse Hoogsteen, or Wobble interactions can be formed by PNA/DNA 
duplexes [l 6 |. 

Unlike DNA which has a helical conformation due to its chiral centers, 
standard PNAs, with no chiral carbon, have no defined structural confor¬ 
mation in solution. Adding a chiral center on the PNA backbone forces left 
or right helix [T7]. and increases PNA/DNA hybridization. Adding a Lysine 
improves water solubility [H], however the Serine is less disruptive to hy¬ 
bridization and -OH groups can form hydrogen bonds with water in order 
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to increase solubility. With the help of chiral centers, PNA homoduplexes 
should also be observed, and eventually present secondary helical structures. 
In this study, we explore the biophysical properties of PNAs built with a chi¬ 
ral center. The chiral center chosen here, consists in a L-Serine substitution 
of the 7 -position on the PNA backbone [H!]. see insert in figure 

In order to determine homo-hybridization of PNAs through a structural 
pre-organization, a series of y-PNAs were analyzed by Circular Dichroism 
(CD) and Liquid State Nuclear Magnetic Resonance (NMR). NMR allows 
to spotlight hydrogen bounds formed by nucleic acid base-pairing. Working 
in neutral conditions, specific NMR experiments bring information about 
the imino hydrogens engaged in the Watson-Crick, Hoogsteen, or Wobble 
bounds. Melting transitions (Tn) of the highlighted PNA homoduplexes 
were determined by CD and absorbance experiments. Circular Dichroism is 
finally used to assess the secondary structure displayed by the homoduplexes. 
All the experimental results show unambiguously that PNAs duplexes are 
formed, and that the chiral centers improve their stability, and drive the du¬ 
plexes into a chiral secondary structure, probably organized as a left-handed 
helix. 

Experimental 

Materials 

PNAs were synthesized as previously reported and purified by HPLCdS]. 
Six PNAs (see figure with a schematic view) were chosen according to their 
nucleobases and the presence or absence of 7 modification on the backbone. 
PNA 1 is a small one, with three nucleobases GGT and has a standard back¬ 
bone. PNA 2 and 3 have the same six nucleobases: GCCGGT, differing by the 
presence of the 7 L-Serine on PNA 3 schematically represented by a star. 
PNA 4 and 5 are also similar PNAs, their backbone contain ten nucleobases 
TGCCGGTTCC, and differ only by the presence of the 7 L-Serine on PNA 4. 
PNA 6 is the complementary strand of PNA 4 and 5, and presents also chiral 
centers. 

First, PNAs were dissolved in Milli-Q water, and diluted into phosphate 
buffer (pH= 6 . 8 ). 10% of deuterated water was added and 10% of deuterated 
DMSO was also added to PNA 5 solution to insure solubility. Before study, 
all solutions were exposed to several annealing cycles by using a solid bath, 
and raising the temperature from room temperature to 95 °C and slowly 
cooled back to room temperature. 
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Figure 1: Schematic view of the 6 PNAs analyzed. Insert shows PNA backbone and the 
7 position of L-Serine. Stars represent chiral centers on the PNA backbone. PNA 2 and 3 
have the same nucleobases sequence but different backbones: 3 is a q-PNA with 3 chiral 
centers. PNA 4 and 5 have the same nucleobase sequence. PNA 6 is the complementary 
PNA of 4, and PNA 7 is the 4-6 duplex. Note that only PNA 5 has a Biotin and not a 
Glutamine at the end of the backbone chain. PNA 1, 2 and 5 have no chiral centers. 

NMR spectroscopy 

PNAs were analyzed in 3 mm tubes on a 700MHz Bruker spectrometer 
equipped with a Z-grad triple resonance cryoprobe. PNAs were prepared 
as above, with 10% D 2 O added for lock. Watergate or Excitation Sculpting 
water suppression experiments were not used because of the concern that the 
exchangeable imino protons would be suppressed too. The ID suppression 
experiments used here is the Jump & Return (JR) excitation sequence pi], 
followed by a W5 water removal sequence [21]. The JR sequence suppresses 
the water signal by not exciting the water spin, in consequence the water 
magnetization is unperturbed and exchangeable protons like the imino pro¬ 
tons appear fully in the ID NMR spectrum. The W5 sequence cleans 
the spectrum from residual water signal, with no perturbation on the other 
signals. 

NMR experiments were run at a various temperatures, between 288 K and 
318 K, as noted in the figure captions. In order to investigate the structural 
conformation and to assign spectra, JR-NOESY, HSQC, TOSCY and DOSY 
spectra were recorded on each compound. Spectra not presented in the text 
can be found in the Supplementary Materials. PNA concentrations range 
from 0.88 mM to 1.71 mM depending on the experiments and the studied 
PNA. 
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Optical Measure 
Absorbanee. 

PNA 3, 4 and 5 were analyzed at 60 /iM (PNA 3) and 50 fjM. (PNA 4 - 
5) by Absorbance on a Jasco J-815 spectropolarimeter in 1 mm path length 
quartz UV cells. 

Tm determinations were performed by absorbance measurement as a func¬ 
tion of the temperature. From an initial absorbance spectrum, five wave¬ 
lengths were chosen and followed during the temperature variation. Ab¬ 
sorbance curves were recorded at several positions between 200 nm and 350 nm 
corresponding to the highest absorbance, or the lowest absorbance (blank 
signal). Absorbance was monitored during heating, the temperature varying 
from 308 K to 368 K with a 1 K/mn variation rate, and a cooling experiment 
performed at the same rate. The were determined by a fitting proce¬ 
dure of the various absorbance and CD curves obtained, using a python 
program (see Supplementary Materials). The program is freely available at 
https://github.com/delsuc/Melting-curve-analysis 

Cireular Dichroism. 

PNA 3, 4, 6, and 7 were analyzed at 60 fjM (PNA 3) and 50 pM (PNA 
4-6-7 ) by Circular Dichroism on the same Jasco equipment. 

Samples were heated from 308 K to 363 K, and cooled down from 363 K 
to 308 K with the same temperature slope. CD spectra were first measured 
at 308 K before the temperature variation, then at 363 K, and finally again 
at the low 308 K after the cooling procedure. 

Melting transitions were compared to prediction obtained for RNA pre¬ 
senting the same sequence, using the DINAMelt Web Server (Di-Nucleic 
Acid hybridization and Melting prediction). |22] This tool predicts nucleic 
acid pairing and melting transitions for DNA or RNA strands. It is found 
at http;//mfold.rna.albany.edu/?q=dinamelt. Schematic view of PNA 
homo-hybridization proposed by the program are shown in Supplementary 
Materials (Figure S14). 

Results and Discussion 

NMR solution analysis of PNAs 

The various PNA constructs were first studied by NMR. The NMR spec¬ 
trum of PNA 3 is shown in figure Because of the JR sequence used, the 
spectrum presents two domains, a positive one and a negative one, on each 
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Figure 2: ID JR NMR spectra of PNA 3. NS=196, [3]=l.lmM, [Phosphate]=25mM, 
T=298K, 700MHz. Signals at 12.9-13 ppm correspond to imino protons involved in G-C 
base pairing, and signals at 12 ppm correspond to imino proton from T-G base pairing. 


side of the location of the water signal. The use of the JR sequence allows 
the observation of the exchangeable protons which would be lost using a 
standard ID excitation. 

By homology with the NMR signals of DNA duplexes, the 13 ppm res¬ 
onance signals are assigned to the imino proton of G-C Watson-Crick base 
pairs. Signals at 12 ppm are assigned to the imino proton of a Wobble base 
pair between a G nucleobase and a T nucleobase. From the NMR signals 
alone, it is not possible to determine whether the spectrum corresponds to a 
duplex formation or to a folding of the PNA on itself in a hair-pin conforma¬ 
tion. The DOSY experiment performed on this molecule has been inconclu¬ 
sive (see Supplementary Data S3). However, the presence of two G-C base 
pairing seems incompatible with a hair-pin geometry, and probably indicates 
that PNA 3 is in a homoduplex form in solution. 

The comparison of the integrals of the imino and amino protons is difh- 
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Figure 3: ID proton JR NMR spectra of 7 -PNA 2 (left) and PNA 3 (right) at vari¬ 
ous temperatures, zoomed in the 11.4-13.4ppm area. [ 2 ]= 0 . 88 mM, [ 3 ]=l.lmM, [Phos¬ 
phate] =25mM, Signals at 12.9-13 ppm correspond to G-C base pairing, and signals at 
12 ppm correspond to T-G base pairing. 

cult because the excitation sequence used here has a non-uniform excitation 
prohle over the spectrum, however it can be estimated from the respective 
integrals that there is a default of imino signals in the spectrum. As imino 
protons are specific to stable base-pairs, this is indicative that some amount 
of the PNA material is still in a single strand form in solution. 

NMR spectra of PNA 2 and 3 for a temperature variation are shown in 
figure The figure presents a zoom of the imino protons area, and spectra 
are recorded for temperatures ranging from 288 K to 318 K. For PNA 3, NMR 
spectra at various temperatures do not present any differences, in contrast 
to PNA 2 for which differences are seen for the T-G imino proton signal 
with a decrease of the signal at high temperature. For PNA 2 the T-G 
signal disappears completely at 318 K, unlike the G-C signals, which are still 
present. This indicates a stronger binding for the G-C base pair compared 
to the T-G base pair. As PNA 2 and PNA 3 contain the same nucleobase 
sequence, this indicates that the presence of the L-Serine in the 7 position 
strengthen the Wobble T-G pair and increases the overall stability of the 
PNA homoduplex. The same temperature variation experiment has been 
made for PNA 4 and 5, (see figures S 12 and S13 Supplementary Material) 
and lead to the same conclusion. 

The NMR instrument cannot go higher in temperature due to the limit on 
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the equipment, and the melting transition is not visible in this experiment. 
Optical approach were thus utilized to determine precisely the melting tran¬ 
sition of the different PNAs. 

Secondary structural analysis 

The melting-transition temperature Tm of the various q-PNA studied here 
have been determined. Melting temperature were determined by tempera¬ 
ture variation, monitoring both the absorbance and the ellipticity at several 
different wavelengths (see Supplementary Materials S17-S35). The experi¬ 
mental values are presented in the table They are compared to theoretical 
Tm values that nucleic acid strands with the same sequence would present, as 
computed using the DINAMelt Web Server Values computed for RNA 
were considered, because the enhanced flexibility of the RNA backbone better 
matches the flexibility of the q-PNA structure. 

Except for PNA 1 which do not present any duplex formation, the differ¬ 
ent PNAs studied here present strong base pairing, with rather high melting 
temperatures. The addition of a chiral center, in the form of a L-Serine 
substitution in the 7 position on the PNA backbone has a stabilizing effect. 
This substitution has no net effect in the case of PNAs 2 and 3 but raises 
by 14 degrees the longer PNAs 5 and 4. The duplex 7, with 10 potential 
base-pairs and a chiral backbone for both strands, presents an high stability, 
with a transition temperature well above 100 °C. Such an extreme stability 
implies that we cannot be sure of the complete pairing of the strands in 
this sample, as the annealing cycles, performed in water, cannot completely 
disrupt the duplex in order to allow a complete sampling of all the possible 
conformations. 

For all the PNAs tested here, the experimental melting temperatures are 
quite higher than the computed values of their RNA counterpart. This is 
probably due to the absence of the strong Coulombic repulsion observed in 
RNA because of the charged backbone. It should be noted that the differences 
between the experimental and computed values present important differences 
for the considered PNAs. In the case of PNAs 2 and 3, with four Watson- 
Crick G-C and 2 Wooble T-G base pairs, this difference is 30 to 35 °G. PNAs 
4 and 5 have the same potential base-paring, located on a longer strand. 
They present differences on the same order of magnitude, except may be for 
PNA 5 which lacks chiral centers, and for which the entropic effect of the 
long non paired strand may have some destabilization effect. PNA 6 on the 
other hand is puzzling. While it is expected to form only 4 Watson-Grick 
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Exp Tm 

computed Tm 

Diff 

PNA 2 

341 K ± 2.0 

308.8 K 

32.2 

PNA 3 

341 K ± 1.8 

308.8 K 

32.2 

PNA 4 

359 K ± 3.4 

323.8 K 

35.2 

PNA 5 

345 K ± 1.9 

323.8 K 

21.2 

PNA 6 

357 K ± 2.8 

298.4 K 

58.6 

PNA 7 

> 383 K 

342.0 K 

> 40 


Table 1: Experimental and computed melting transition temperatures. 


G-C base pairs, it presents a melting transition equivalent to PNA 4, 58 °C 
above its theoretical RNA counterpart. This cannot be explained without 
the formation of additional non-standard base pair, such as the C-»-A or 
A-»-A base pairs, as observed in some RNA secondary structures [23]. These 
non-standard pairs are not isosteric to the standard Watson-Crick ones, but 
their formation might be possible here thanks to the greater flexibility on the 
PNA backbone. 

Full circular dichroism spectra of PNA 3 and 4 recorded at 308 K and 
368 K are presented in figure The ellipticity signals are very different at 
low and high temperatures. At 308 K (red curve) the complex signal observed 
for both PNAs is characteristic of a secondary structure existing in solution. 
The dichroic signal varies according to the chirality of the molecule and of 
the secondary structure of the molecule. This shows a pronounced Cotton 
effects, characteristic of a helix as show by Dragulescu-Andrasi et alUT]. 
This secondary structure is disrupted at high temperature (green curve), as 
at 368 K the ellipticity signal presents much less structure and remains mostly 
flat, probably dominated only by the chirality of the backbone C.y carbons. 
This absence of secondary structure at high temperature is characteristic of 
the melting of the duplex, as was already indicated by the hyperchromicity 
observed at high temperatures (see Supplementary Material figure SI7). The 
curves recorded at 308 K before and after (blue curve) the high temperature 
denaturation are fully superimposed, indicating that the secondary structure 
of PNA in solution is fully reversible. All these results are indicative of the 
formation of a chiral supramolecular organization, probably in the form of a 
helical secondary structure of the double strand. 

Comparison of the overall shape of the CD spectra with the literature [21], 
and considering an antiparallel duplex, the maxima around 262 and 217 nm. 
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Figure 4: Ellipticity curves of PNA 3 (left) and PNA 4 (right) at three temperatures. Red 
and blue curves are ellipticity signal at 308 K before and after denaturation respectively, 
green curve is ellipticity measured at 368 K. 



Figure 5: Ellipticity curves of PNA 4 (blue), PNA 6 (green) and PNA 7 (red) as a function 
of the wavelength. PNA 7 is duplex of PNA 4 and PNA 6. 


and the minima around 277, 238, and 200 nm indicate that the formed duplex 
is probably a left handed helix. 

In the same manner, ellipticity curves for PNA 4, 6 and 7 are presented 
in figure All three structures display similar curves, indicating similar 
supramolecular organization. However upon detailed analysis, we can observe 
that the maxima of the CD spectra occur at slightly different wavelengths, 
characteristic of slightly different structural organizations. 

Conclusion 

In conclusion, we have shown that PNA can adopt homo- and hetero¬ 
duplex conformations in solution through the formation of standard Watson- 
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Crick as well as non-standard base pairing. These duplexes are extremely 
stable and present melting transitions at temperatures higher than their nu¬ 
cleic acid homologues. An extreme melting transition well above 100 °C was 
observed for a 10 bases complementary duplex, more than 40 °C higher than 
predicted for an RNA equivalent duplex. 

PNAs duplexes with a chiral backbone present a marked chiral secondary 
structure, indicating that they are organized in supramolecular helices. Cir¬ 
cular Dichroism spectra indicate a common folding pattern for all studied 
structures, with nevertheless small differences depending on the details of 
the nucleobase sequence. 
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Abstract 

This documents presents additional spectra obtained on the different PNA 
molecules. 
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Figure SI: ID NMR spectrum of PNA 1 at 700MHz with solvent presaturation. 
[l]=1.56mM, [Phosphate] =25mM, T=298K. 

NMR spectrum of PNA 1 shows more aromatics protons than expected, 
it seems to have several methyls groups. This sample is not completely clean. 
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Figure S2: 2D ^^C-HSQC NMR spectrum of PNA 1 at 700MHz. [l]=1.56mM, 

[Phosphate] =25mM, T=298K. 
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Figure S3: 2D DOSY NMR spectrum of PNA 1 at 700MHz with excitation sculpting. 
[l]=1.56mM, [Phosphate] =25mM, T=298K. 


On the HSQC spectrum of PNA 1 (S2), we highlight the presence of 
two or three methyl groups. The PNA 1 is GOT made of, only one methyl 
of thymine should be present. Guanine groups are present on the HSQC 
experiment. As seen on the aliphatique area (3-5ppm), a mix of PNA could 
be present. The DOSY experiment shown in S3 doesn’t confirm the presence 
of other PNA. It is not exclude to have other 3 nucleobases PNA with a 
similar molecular weight and apparent diffusion coefficient. 
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Figure S4: ID Jump & Return NMR spectrum of PNA 2 at 700MHz. [2]=880/uM, 
[Phosphate] =25mM, T=298K. The DINAMelt prediction for homodimerization is drawn. 


PNA 2 should form homodimer as drawn in the insert of the pictnre S4. 
G-T and G-C base pairing are apparent on the Jump & Return NMR spec¬ 
trum, indicating a self base pairing. Intensities of peak are non proportional 
to the proton nnmber in the Jump & Return experiment as soon as they are 
far from the center of the window. Thns, the G-C imino protons peak appears 
weaker than the G-T imino proton peak which have twice a less number of 
protons. 
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Figure S5: 2D ^ H Jump & Return NOESY NMR spectrum of PNA 2 at 700MHz. 
[2]=880/uM, [Phosphate] =25mM, T=313K, mixing time= 250msec. 


The NOESYjr NMR spectrum of PNA 2 does not highlight presence of 
cross peak between G-C imino proton and water. 
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Lgure S6: Partial assignment of PNA 3 at 700MHz. [3]=l.lmM, [Phosphate] =25mM, 
=298K. 





Figure S7: 2D ^^C-HSQC NMR spectrum of PNA 3 at 700MHz. [3]=l.lmM, 

[Phosphate] =25mM, T=303K. Acquisition time= 1 day. 
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Figure S8: Zoom of 2D H Jump & Return NOESY NMR spectrum of PNA 3 at 

700MHz. [3]=l.lmM, [Phosphate] =25mM, T=303K, mixing time= 250m5ec, acquisition 
time= 1 day and 2 hours. 
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PNA3 (ID IH) 














Figure S9: Zoom of 2D H TOCSY NMR spectrum of PNA 3 at 700MHz. 

[3]=l.lmM, [Phosphate] =25mM, T=303K, mixing time= 90msec, acquisition time= 19 
hours. 
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Temperature variation of PNA 4 shows at least two G-C and two or three 
G-T imino protons, The G-T base pairing appears strong as G-T imimo 
protons signal remain visible at 318K. 
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Figure SIO: ID Jump & Return NMR spectrum of PNA 4 at 700MHz. [4]=1.16mM, 
[Phosphate] =25mM, variable temperature. 
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In comparison the to PNA 4, the Biotin ended PNA5 shows one G-C 
base pairing and one large signal for a G-T base pairing. The G-T base 
pairing appears less stable as the G-T imino proton signal decrease at high 
temperature. This imino proton signal was still evident on PNA 4 at 318K, 
whereas this signal vanish at the same temperature in PNA 5, meaning the 
biotin impact on the folding of these PNA. 
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Figure S12: ID Jump & Return NMR spectrum of PNA 5 at 500MHz. [5]=1.71mM, 
[Phosphate] =25mM, variable temperature. 
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Figure S13: 2D H Jump & Return NOESY NMR spectrum of PNA 5 at 700MHz. 

[5]=1.71mM, [Phosphate] =25mM, T=313K, mixing time= 4:00msec. 


For the first time, the NOESY spectrum of PNA 5 shows correlation 
between one imino proton and protons from the backbone chain. Sadly, the 
complete assignement of this PNA remain impossible on this sample. 
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Figure S14: DINAMelt prediction of homodimerization of PNA 5 and 4 (right) and 6 
(left). 
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Melting_analysis 


August 12, 2015 


1 Duplex formation and secondary structure of 7 -PNA observed 
by NMR and CD Supplementary Material 

2 Content 

This file present the program used to determine the melting temperatures Tm of PNA molecules, from 
absorbance and CD data extracted acquired on Jasco spectrometer. 

The document is in two parts 

• the fitting program itself, 

• the results obtained when applying the program to experimental data provided in another file-set. 
This Program are provided under the Licence CeCILL 2.1 

This program HAS NOT been tested intensively, it is believed to do what it is supposed to do, Howver, 
you are welcome to check it on your own data. 

• Author : M-A Delsuc (madelsuc@unistra.fr) 

• Date : August 2015 

• Version : 2.0 (Multiexp and error bars was added from version 1.0, also some modification have been 
made to make the prgm easier to use ) 

The program and the datasets are freely available at https://github.com/delsuc/Melting-curve-analysis 


3 The program 

This program is a general fitting program, meant to analyze melting curves, either from Absorption or CD 
data. It has been used here to analyze the PNA data presented in the associated work, but can certainly be 
used on any other data. 

It is in three parts : 

• reading the Jasco files 

• This part in independent from the other part, and can be extended to adapt your own set-up, the main 
entry point is loadO . 

• JASCO CSV and txt files are programmed 

• the FrenchCoding flag allows to access the files on a french version of Windows, 

• defining a theoretical melting curve, used for fitting 

• this is supposed to be general to any copperative melting mechanism 

• The program itself that read files, fit the melting curve and display the results. 

The program comes in two versions: 

• MonoexpO permits to fit one melting experiment, measured at one given wavelength 

• MultiexpO permits to fit several melting experiments, measured at several wavelengths 
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In [1] : #fiTst, set-up the python scene 
import sys 
import codecs 
import numpy as np 
import matplotlib as mpl 
import matplotlib.pyplot as pit 
from scipy.optimize import curve_fit 
Debug = False 

# these seaborn is optional, used for a nice setup 
try: 

import seaborn as sns 
current_palette = sns.color_palette() 
except ImportError: 
pass 

yomatplotlib inline 

mpl.rc( "figure" , figsize=(12, 6)) 

from IPython.display import HTML, display_html 
TheFig = 15 
def Fig(caption=None): 
global TheFig 

St = "<b>Figure syod</B>"yoTheFig 
TheFig += 1 
if caption: 

St += " : <i>yoS</i>"yoCaption 

display_html(HTML(st)) 

3.1 utilities to read the Jasco files 

Remark Our JASCO system is running on a french version of Windows. Files are coded with the latin_l 
coding, depending on your set-up, you may need to change this parameter if you experiment coding errors 
while reading the file. 

On the french system, decimal numbers are coded with a comma for decimal : 33,2 instead of 33.2 
Finally the . csv files on the french system are coded with a semi colon to separate the field, so a line in 
a CSV file appears as : 

1; 33,2; 64,5 

instead of 

1, 33.2, 64.5 

This is taking care of with the FrenchCoding variable set to True 

In [2] : # adapt FileCoding to the ASCII coding used in the JASCO text files 

# exemples are 

# ”cpl252” (Western Europe) ”latin_l” (France and other) ”shift_jis” (Japanese) ”utf_8” 

# check https://docs.python.org/2. 7/library/codecs.html^standard-encodings 
^FileCoding = ’•utf_8” 

FileCoding = "latin_l" 

# set to FrenchCoding to False if you the cvs and txt files are not using the french coding 

# C for decimal point and ; for field separator ) 


2 


#FrenchCoding = False 
FrenchCoding = True 

#Tead files 

def loadcsv(f ich): 

// n n 

reads a csv jasco file 

returns data, meta 

where data as a numpy array 

meta as a dictionary key ’.value 

It It It 

data = [] 
meta = {} 
xydata = False 

for lin in codecs.open(fich, "rb", FileCoding): # the codec may he different on your set-u 

if FrenchCoding: 

lin = lin.replaceC ^ \ ^. O 
linspl = lin. stripO . split ( ^ ; O 
else : 

linspl = lin.strip 0.split( ^, O 
if linspl == [^^]: # marks end of data 

xydata = False 

if xydata: # reading data 

flin = [ float (i) for i in linspl] 
data.append(flin) 

elif len(linspl) >1 : ^ ou pas un champs a recuperer 

meta[linspl[0]] = ";" .join(linspl[1:]) 
if linspl == ["XYDATA"] : # marks beginning of data 

xydata = True 

return np.array(data), meta 
def loadtxt (fich): 

II It It 

reads a text jasco file 

returns data, meta 

where data as a numpy array 

meta as a dictionary key ’.value 

II II II 

data = [] 
meta = {} 
xydata = False 

for lin in codecs.open(fich, "rb", FileCoding): # codec may vary 
if FrenchCoding: 

lin = lin.replaceC \ ^^. O # This is for french coding on windows ! 
lin = lin.stripO 
linspl = lin.split( ^\t O 

if linspl == [^^]: # marks end of data 

xydata = False 

if xydata: # read data 

flin = [ float (i) for i in linspl] 
data.append(flin) 

elif len(linspl) >1 : # if meta 

meta[linspl[0]] = "\t " .join(linspl[1:]) 
if linspl == ["XYDATA"] : # marks beginning of data 
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xydata = True 

return np.array(data), meta 
def load(fich) : 

"read jasco file, txt or csv based on extension .csv" 
try: 

if f ich. endswithC . csv") : 

return loadcsv(fich) 
else : 

return loadtxt(fich) 
except: # generic error message. 

if FrenchCoding: 

decimal = "," 
else : 

decimal = "." 
print """ 

The data file could not be read 
please check 

- file name : 

- file coding : Xs 

- decimal coding : 

"""yoCfich, FileCoding, decimal) 
sys.exit 0 

def decoupe (d): 

"extract columns from data array" 

Ido = d[:,0] 
cd = d [: , 1] 
volt = d [:,2] 
absorb = d [:,3] 

return (Ido, cd, volt, absorb) 

def readvalues (fich, tampon=None, lincor=False) : 

// // // 

read file and return cd and ahs 

if tampon is a file, will corrected for buffer 

if lincor is True, a linear component is removed 

returns lambda_array, absorb_array, cd_array, remark_text 

n n n 

d,meta = load(fich) 

# print meta[’Comment’] 

(Ido, cd, volt, absorb) = decoupe(d) 
if tampon: 

dtamp,metatamp = load(tampon) 

(a, b ,c , absorb_tamp) = decoupe(dtamp) 
absorb = absorb-absorb_tamp 
rem = " - corrected" 
else : 

rem = "" 

if lincor: ^eventually correct for linear component evaluated on first points 

zone = int (0. l*len(ldo) ) 

fit = np.polyfit(ldo[0:zone], absorb[0:zone], 1) 
absorb = absorb-np.polyval(fit, Ido) 
return Ido, absorb, cd, rem 
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def plotcd(fich, tampon=None , lincor=False , label=None) : 

"draws cd from file:fich with label, with eventual corrections" 

Ido, absorb, cd, rem = readvalues(fich, tampon=tampon, lincor=lincor) 
if label is None: 

label = fich+rem 
pit.plot(Ido, cd, label=label) 
pit.legend0 
pit.titleCfich) 

def plotabs (fich, tampon=None, lincor=False , label=None) : 

"draws absorbance from file:fich with label, with eventual corrections" 

Ido, absorb, cd, rem = readvalues(fich, tampon=tampon, lincor=lincor) 
if label is None: 

label = fich+rem 

pit.plot(Ido, absorb, label=label) 
pit.legend0 
pit.title(fich) 

3.2 melting curve definition, 

used for fitting, implements the hypoerbolic-tangent model, see for instance : M. YA. AZBEL, “DNA 
sequencing and melting curve” P.N.A.S. 76 (1) pp.101-105 (1979) 

In [3] : #r_/7? fitter 

def Tmfuncd, lo, le, Tm, beta): 

// // It 

draw a theoretical fusion curve 
T: temperature range 
lo le : value at low and high temp 
Tm: fusion temperature 

beta: cooperativity parameter, proportional to the inverse of the temperature range over wh 

n n n 

if Debug: 

print "T",T,"\n args",(Io, le, Tm, beta) 
tred = beta*(T-Tm) 

return lo + 0.5*(le-Io)* (1+np .tanh(tred)) 

#and plot theoretical curve 
X = np.Iinspace(0,100,100) 

plt.plot(x, Tmfunc(x, 5,10,60,0.1), label=r "$\beta=0. 1$") 
pit.plot(x, Tmfunc(x, 5,10,60,0.07), label=r "$\beta=0.7$" ) 
pit.plot(x, Tmfunc(x, 5,10,60,0.05), label=r"$\beta=0 .05$" ) 
pit .xlabeK "Temperature ($''\circ C$)") 
plt.ylabelC'A.U.") 
pit.Iegend(loc=0) 

Fig(r"Theoretical melting curves for a $T_m$ of $60"\circ C$ with varying $\beta$") 

Figure S15 : Theoretical melting curves for a $Tjii$ of $60"\circ C$ with varying $\beta$ 
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3.3 the fitting program 

The class MonoexpO defines an object which loads the various data of a temperature run followed by 
CD/Absorbance at one wavelength. 

To use, initialize the various attributes (parameters) (check for documentation), attributes are : 

• BASE : the base directory which holds the files 

• fnom : the generic name for the file 

• freq : the frequency (wavelength) on wich to apply the measure 

• cell : the cell to follow 

• tampon : None, or the cell of the buffer 

• meth : meth is either “ABS” or “CD” 

• nom : used for label in plots 

• extension : “.txt” or “.csv” 

then use either : 

• plotO : plots the curve 

• analyze 0 : plot and fit melting curve, determines and p 

In [4] : ##### Create the •programs that compute and display results 

class Monoexp(object) : 

// // // 

This class defines an object which loads the various data which describe one experiment 
of a temperature run followed at one wavelength by CD/Absorbance. 

To use, 

- initialize the various parameters (check ^ _ init _ () ^ function for detained documentati 

then use either : 

- plotO : plots the curve 

- analyzeO : plot and fit melting curve 
// // // 
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def _init_ (self) : 

// // // 


creates the object 

The file name will he built from these parameters with self.full () 

It n It 


self. BASE = 

self.fnom = "filename-" 

self .freq = 217 

self. cell = 2 

self. tampon = None 

self.meth = "abs" 

self.nom = None 

self .extension = ".txt" 


# the base directory which holds the files 

# the generic name for the file 

# the frequency to measure 

# the cell to follow 

# None, or the cell of the buffer 

# meth is either "ABS" or "CD" 

# used for label in plots 

# ".txt" or ".CSV" 


def Is(self): 

"utility to print all possible files, using the name preset" 
import glob 

fn = "yos/yoS*yoS "% (self . BASE, self.fnom, self . extension) 

print "List of files :" 

print "\n".join( glob.glob(fn) ) 

©property 
def full (self): 

"the filename" 

return "yos/yosyodnm-Cell ydyoS "% (self . BASE, self.fnom, self, f req, self, cell, self . extensio: 
©property 
def tfull(self) : 

"the filename of the buffer" 
if self. tampon is not None: 

return "yoS/yoSyodnm-Cell ydyoS "^(s elf . BASE, self.fnom, self, f req, self, tampon, self, ex' 
else : 

return None 
©property 
def title(self) : 

"title for plot" 

if self .meth.lower0 == ^cd^: 

return ^CD over Temperature^ 
elif self .meth.lower0 == ^abs^: 

return ^Absorbance over Temperature^ 

©property 

def label(self) : 

"label for plot" 
if self.nom is None: 

self.nom = self. cell 

return ^PNA Xs - $\lambda=yod\,nm$^ %( self.nom, self. freq) 

def fit (self, Tm=60.0, beta=0.1, Io=None, Ie=None) : 

// // // 

fit and plot result, initial estimate can be given to help the fit 

uses all current attributes (see _ init _ ) to determine which dataset is to he fitted 

results are stored in the self?Results dictionary 
// // // 

global Debug 

Ido, absorb, cd, rem = readvalues (self .full, tampon=self .tfull) 
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self. Ido = Ido 

if self .meth.lower0 == "abs": 
y = absorb 

elif self .meth.lower0 == "cd": 
y = cd 

# set-up initial fit values 
if lo is None: 

lo = y[0] 
if le is None: 
le = y[-l] 

try: 

popt, pcov = curve_fit(Tmfunc, Ido, y, pO=[Io, le, Tm, beta] ) # fit 
except RuntimeError : 

print self. label, "curve could not be fitted" 
self. Results = None 
else : 

perr = 2*np.sqrt(np.diag(pcov)) 
self. Results = { 

^lo ^ :popt [0], 

^le ^ :popt [1], 

:popt [2], 

^beta^ :popt [3], 

^Tm_errorbar^ :perr[2], 

^beta_error^ :perr[3] }# store monofit 
self.residu = y-Tmfunc(Ido, *popt) 
self. norm = np.sqrt(np. sum(self .residu**2)) 
if Debug: print pcov 

print "“/oS Tm: “/.IfC +/- %.2± beta: %.'1± chi2:yo.3f "70 (self .label, popt [2] , perr[2] 
def showfit (self ) : 

"plot the result of the fit" 
try: 

R = self. Results 
except : 

raise ExceptionC'f it is not performed yet") 
if R is not None: # None indicates fit not converged 

pit.plot (self .Ido, Tmfunc (self .Ido, R["Io"] , R["Ie"] , R["Tm"] , R["beta"]), ^k— O # 
mn = min( R["Io"] , R["Ie"] ) 
mx = max( R["Io"] , R["Ie"] ) 

if R["Tm"] < 100.0 and R["Tm"] >0.0: # do not plot out-of-range results 

plt.plot([R["Tm"] ,R["Tm"]] , [mn,mx] , 

def plot (self): 

"plot the experimental curve" 
if self .meth.lower0 == "abs": 
self .abs 0 
plt.ylabelC'A.U.") 
elif self .meth.lower0 == "cd": 
self .cd() 

pit .ylabelC'Ellipticity - millideg") 
pit.Iegend(loc=0) 
pit .title(self .title) 
pit .xlabeK "Temperature ($''\circ C$)") 
def analyze (self , Tm=60.0, beta=0.1, Io=None, Ie=None) : 

"fit (using initial guess if provided) and plot" 
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self .plot 0 

self .fit(Io=Io, Ie=Ie, Tm=Tm, beta=beta) 
self .showfit 0 
def abs(self): 

"read absorbance values" 

Ido, absorb, cd, rem = readvalues (self .full, tampon=self .tfull) 
pit.plot(Ido, absorb, label=self .label) 
def cd(self ): 

"read CD values" 

Ido, absorb, cd, rem = readvalues (self .full, tampon=self .tfull) 
pit.plot(Ido, cd, label=self .label) 

3.4 second method 

The class MultiexpO is very close to Monoeexp, the only difference, is that it allows to fit several measures 
at different wavelength, to the same Tm value 

It is to initialize in the same manner than Monoexp, It has an additional attribute : 

• multifreq : a list of the wavelengths to use in the fitting/plotting operations 
then use either : 

• plotO : plots the curve 

• analyze 0 : plot and fit all the curves named by multifreq to a same 

In [5]: def TmMfunc(T, *args): 

// // // 

Equivalent to Tmfunc, bur for several wavelengths, using the same Tm and beta 
everything is flattened into ID arrays (parameters, and return values) 

T: temperature range 
N, lo, le, Tm, beta = *args 
Tm: fusion temperature 
beta: cooperativity parameter 
N: number of curve to analyze, 
then for each curve 

To le : value at low and high temp 

It // // 

global Ncurves, Debug # Ncurves contains the number of curves to fit 
if Debug: 

print "T" ,T, " \n args" ,args 
N = Ncurves 
lo = args[0:N] 
le = args[N:2*N] 

Tm = args[-2] 
beta = args[-l] 
r = [] 

for i in range(len(Io) ): 
tred = beta*(T[i]-Tm) 

r.append(lo[i] + 0.5*(le [i]-lo[i])*(1+np.tanh(tred))) 
return np.array(r).flatten() 
class Multiexp (Monoexp): 

def _init_ (self) : 

super (Multiexp, self). _init_() 

self .multifreq = [] 
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def fit (self, Tm=60.0, beta=0.1, Io=None, Ie=None) : 

// // // 

fit and plot result, initial estimate can he given to help the fit 
values are given for initial estimate 

uses all current attributes (see _ init _ ) to determine which dataset is to he fitted 

results are stored in the self.Results dictionary 
// // // 

global Ncurves # Ncurves contains the number of curves to fit 

y = [] 

self. Ido = [] 

for self.freq in self .multifreq: 

Ido, absorb, cd, rem = readvalues (self .full, tampon=self .tfull) 

self .Ido.append(Ido) 

if self .meth.lower0 == "abs”: 

y.append(absorb) 
elif self .meth.lower0 == "cd": 
y.append(cd) 

# set-up initial fit values 
if lo is None: 

lo = [i [0] for i in y] 
if le is None: 

le = [i[-l] for i in y] 
y = np.array(y).flatten0 
self. Ido = np.array (self .Ido) 

Ncurves = len(Io) 

pO = lo + le + [Tm] + [beta] # all parameters are concatenated 
r = (TmMfunc (self .Ido, *p0 ) - y) 
if Debug: 

print self .Ido.shape, y.shape, len (pO) 
print TmMfunc (self .Ido, *p0 ).shape 
print pO 
print r 

try: 

popt, pcov = curve_fit(TmMfunc, self. Ido, y, p0=p0, maxfev = 6000) 
except RuntimeError : 
popt = pO 
perr = None 
else : 

print self .fitlabel, "curve could not be fitted" 
if Debug: print popt 

self.residu = y-TmMfunc (self .Ido, *popt) 
self. norm = np.sqrt(np. sum(self .residu**2)) 
perr = 2*np.sqrt(np.diag(pcov)) 

print "yoS Tm: ydfC +/- yo.2f beta: yo.2f chi2 :%.3f "“/(self .f itlabel, popt [-2] , pe; 
if Debug: 

print popt 

self.residu = y-TmMfunc (self .Ido, *popt) 
self .norm = np.sqrt(np. sum(self .residu**2)) 
if perr is not None: 
self. Results = { 

^N^ :Ncurves, 

^lo ^ :popt [0:Ncurves], 
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^le ^ :popt [Ncurves:2*Ncurves], 

:popt[-2], 

^beta^ :popt[-1], 

^Tm_errorbar^ :perr[-2], 

^beta_error^ :perr[-1] }# store monofit 

else : 

self. Results = None 

def showfit (self ): 

"plot the result of the fit" 
try: 

R = self. Results 
except : 

raise ExceptionC'fit is not performed yet") 
if self. Results is not None: 
for i in range(R[^N^] ): 

pit.plot (self .Ido [i], Tmfunc (self .Ido [i], R["Io"] [i], R["Ie"] [i], R["Tm"] , R["b( 
mn = min( min(R["Io"]), min(R["Ie"]) ) 
mx = max( max(R["Io"]), max(R["Ie"]) ) 

if R["Tm"] < 100.0 and R["Tm"] >0.0: # do not plot out-of-range results 

plt.plot([R["Tm"] ,R["Tm"]] , [mn,mx] , 
def plot (self): 

for self.freq in self .multifreq: 
super (Multiexp, self).plot() 

©property 

def f itlabeKself ) : 

"label for fit" 
if self.nom is None: 

self.nom = self. cell 
return ^PNA yoS^yo( self.nom) 


4 Experimental Results 

A directory called PNA_DATA/ should be available along with this program. It contains a set of files holding 
the data acquired on the various PNA studied. 

The experiments were acquired on by Absorbance and CD on a Jasco J-815 spectropolarimeter using 
1mm path length 110 quartz UV cells. Each run was performed first by heating, the temperature varying 
from 308 K to 368 K with a 1 K/mn variation rate, and then by a cooling experiment performed at the same 
rate. 

Each file contains one measure on one PNA, at one wavelength, on a temperature series. The files are 
distributed among several folders corresponding to different experimental runs. The name of the file codes 
for the sample, the wavelength and the temperature variation (going up or down) 

4.1 PNA 2 and 3 

melting followed by absorbance and CD 

• buffer is in 1 

• PNA 2 in 2 

• PNA3 in 3 


4.1.1 CD Spectrum 

Eirst we plot the CD spectrum : 
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In [6]: plotcd(fich="PNA_DATA/Dichro_130627/PNA8xx_T35a_lmm_130627-3.txt", 

tampon=''PNA_DATA/Dichro_130627/PNA8xx_T35a_lmm_130627-l.txt", lincor=False , label="PNA 3 
plotcd(fich="PNA_DATA/Dichro_130627/PNA8xx_AftDenatT35a_lmm_130627-3.txt", 

tampon="PNA_DATA/Dichro_130627/PNA8xx_T35a_lmm_130627-l.txt", lincor=False , label="PNA 3 
Fig(r"CD spectra of PNA 3 at 308K before and after the variation temperature") 

Figure S16 : CD spectra of PNA 3 at 308K before and after the variation temperature 
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4.1.2 Temperature variation - going-up 

1°C / min 

The following wave-lengths were recorded : 200nm 217nm 230nm 260nm 310nm 380nm 
Results are shown for Absorbance signal, and for CD signal for PNA-3 (PNA-2 is non chiral) 
Not all frequencies give fittable results. 

First we can use MonoexpO to fit the best curve. 

In [7] : f = MonoexpO 

f.BASE = "PNA_DATA/Dichro_130627" 
f.fnom = "PNA8xx_Denat_130627-l-" 
f .IsO 

List of files : 

PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-200nm-Cell 1 .txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-200nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-200nm-Cell 3.txt 
PNAJ)ATA/Dichro_130627/PNA8xx_Denat_130627-l-217nm-Cell 1 .txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-217nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-217nm-Cell 3.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-230nm-Cell 1 .txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-230nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-230nm-Cell 3.txt 
PNAJ)ATA/Dichro_130627/PNA8xx_Denat_130627-l-260nm-Cell 1 .txt 
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PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-260iim-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-260iim-Cell 3.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-310iim-Cell 1. txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-310nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-310iim-Cell 3.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-380iim-Cell 1. txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-380iim-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Denat_130627-l-380iim-Cell 3.txt 


4.1.3 PNA-2 Absorbance 

First we fit one wavelength : 200nni 

In [8]: f.tampon = 1 
f.cell = 2 
f.meth = "ABS” 

f.freq = 200 # the 200nm curve give the best results 

f.analyze 0 

Fig( ''PNA-2 Absorbance") 

PNA 2 - $\lambda=200\,nm$ Tm: 66.2C +/- 0.94 beta: 0.06 chi2:0.014 

Figure S17 : PNA-2 Absorbance 



Then we try to fit all the curves at once using Multiexp(), defined in multifreq 

In [9] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130627" 
f.fnom = "PNA8xx_Denat_130627-l-" 
f.tampon = 1 
f.meth = "ABS" 

f.multifreq = (200, 217, 230, 260,310,380) 
f.analyze 0 

Fig( "PNA-2 Absorbance") 
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PNA 2 curve could not be fitted 

PNA 2 Tm: 66.3C +/- 1.30 beta: 0.06 chi2:0.048 

Figure S18 : PNA-2 Absorbance 
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As MultiexpO gives satisfactory results, we will use it in the following, all curves are plotted, but only 
the curve giving interpretable results are fitted. 

Results are accumulated into a global ordered dictionary, for the final conclusion 

In [10]: import collections 

RESULTS = collections.OrderedDict0 
RESULTS ["PNA-2 Abs"] =f.Results 

4.1.4 PNA-3 Absorbance 
In [11]: f.cell = 3 

f.multifreq = (200, 217, 230, 260,310,380) 
f.analyze 0 

RESULTS ["PNA-3 Abs"] = f.Results 
Fig ("PNA-3 Absorbance") 

PNA 2 curve could not be fitted 

PNA 2 Tm: 71.1C +/- 1.77 beta: 0.04 chi2:0.054 

Figure S19 : PNA-3 Absorbance 
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4.1.5 PNA-3 CD 


In [12]: f.meth = "CD" 
f.cell = 3 

f.multifreq = (200, 217, 230, 260, 310, 380) 
f.analyze 0 

RESULTS ["PNA-3 CD"] =f.Results 
Fig("PNA-3 CD") 


PNA 2 curve could not be fitted 

PNA 2 Tm: 68.4C +/- 0.98 beta: 0.06 chi2:1.572 
Figure S20 : PNA-3 CD 
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4.1.6 Temperature variation, going down 

(1 degree / minute) 

following wave-lengthes were recorded : 217nm 260nm 380nm 

In this series, only CD of PNA-3 could be reliably analyzed, in particular because of the lack of measure¬ 
ments at 200nm, wave-length that gave the best results in the previous series. 

Nevertheless, some shift in temperature are observed compared to the previous experiment. 

In [13] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130627" 
f.fnom = "PNA8xx_Cooling_130627-l-" 
f .IsO 

List of files : 

PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-217nm-Cell l.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-217nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-217nm-Cell 3.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-260nm-Cell l.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-260nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-260nm-Cell 3.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-380nm-Cell l.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-380nm-Cell 2.txt 
PNA_DATA/Dichro_130627/PNA8xx_Cooling_130627-l-380nm-Cell 3.txt 

4.1.7 PNA-2 Absorbance 

In [14]: f.meth = "ABS" 
f.tampon = 1 
f.cell = 2 

f.multifreq = (217, 260, 380,) 
f.analyze 0 

Fig ("PNA-2 Absorbance") 

PNA 2 curve could not be fitted 

PNA 2 Tm: 54.3C +/- 36.58 beta: 0.03 chi2:0.030 

Figure S21 : PNA-2 Absorbance 
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4.1.8 PNA-3 Absorbance 

In [15]: f.meth = "ABS" 
f.tampon = 1 
f.cell = 3 

f.multifreq = (217, 260, 380,) 
f.analyze 0 

Fig ("PNA-3 Absorbance") 


PNA 2 curve could not be fitted 

PNA 2 Tm: 73.2C +/- 5.00 beta: 0.11 chi2:0.025 

Figure S22 : PNA-3 Absorbance 
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4.1.9 PNA-3 CD 


In [16]: f.meth = "CD" 
f.tampon = 1 
f.cell = 3 

f.multifreq = (217, 260, 380,) 
f.analyze 0 

RESULTS ["PNA-3 CD 2"] = f.Results # results stored as 2 are for decreasing temperatures 
Fig("PNA-3 CD") 

PNA 2 curve could not be fitted 

PNA 2 Tm: 68.1C +/- 1.82 beta: 0.06 chi2:0.788 

Figure S23 : PNA-3 CD 
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4.2 PNA 4 & 6 

The same protocole was then applied on the two PNA 4 and 6. Here PNA-4 is in cell 3, and PNA-6 is in 

cell 4. 

PNA-7 (see below) is in cell 5. 

4.2.1 going up 

raison the temperature, 1°C / min 

measured wavelengths 218 240 262 280 350 

In [17] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130711" 
f.fnom = "PNA8xx_Denat_130711-l-" 
f .IsO 
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List of files : 

PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 2.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 4.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat.l30711-l-240nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 5.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat.l30711-l-280nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell l.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat.l30711-l-350nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 5.txt 


4.2.2 PNA-4 Absorbance 

In [18]: f.meth = "ABS" 
f.tampon = 1 
f.cell = 3 
f.nom = 4 

f.multifreq = (218, 280, 240, 262, 350) 
f.analyze() 

RESULTS ["PNA-4 Abs"] = f.Results 
Fig ("PNA-4 Absorbance") 

PNA 4 curve could not be fitted 

PNA 4 Tm: 86.4C +/- 4.18 beta: 0.11 chi2:0.059 

Figure S24 : PNA-4 Absorbance 
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4.2.3 PNA-6 Absorbance 


In [19]: f.cell = 4 
f.nom = 6 

f.multifreq = (218, 280, 240, 262, 350) 
f.analyze 0 

RESULTS [''PNA-6 Abs"] = f. Results 
Fig ("PNA-6 Absorbance") 


PNA 6 curve could not be fitted 

PNA 6 Tm: 88.8C +/- 5.51 beta: 0.09 chi2:0.060 

Figure S25 : PNA-6 Absorbance 
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4.2.4 PNA-4 by CD 

In [20]: f.meth = "CD" 
f.cell = 3 
f.nom = 4 

f.multifreq = (218, 280, 240, 262, 350) 
f.analyze 0 

RESULTS ["PNA-4 CD"] = f.Results 
Fig("PNA-4 CD") 

PNA 4 curve could not be fitted 

PNA 4 Tm: 81.7C +/- 2.43 beta: 0.04 chi2:0.750 

Figure S26 : PNA-4 CD 
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4.2.5 PNA-6 by CD 

In [21]: f.cell = 4 
f.nom = 6 

f.multifreq = (218, 280, 240, 262, 350) 
f.analyze 0 

RESULTS ["PNA-6 CD"] = f.Results 
Fig("PNA-6 CD") 

PNA 6 curve could not be fitted 

PNA 6 Tm: 83.OC +/- 4.50 beta: 0.03 chi2:1.108 
Figure S27 : PNA-6 CD 
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4.2.6 going down 

218 240 262 280 350 Abs does not work 

In [22] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130711" 
f.fnom = "PNA8xx_Cooling_130711-l-" 
f .IsO 

List of files : 

PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-218nm-Cell 1 .txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-218nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-218nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-218nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-218nm-Cell 5.txt 
PNA_DATA/Dichro.l30711/PNA8xx_Cooling_130711-l-240nm-Cell 1 .txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-240nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-240nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-240nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-240nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-262nm-Cell 1 .txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Cooling.l30711-l-262nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-262nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-262nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-262nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-280nm-Cell 1 .txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-280nm-Cell 2.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Cooling.l30711-l-280nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-280nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-280nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-350nm-Cell 1 .txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-350nm-Cell 2.txt 
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PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-350nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-350nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Cooling_130711-l-350nm-Cell 5.txt 

In [23]: f.meth = "CD" 
f.tampon = 1 
f.cell = 3 
f.nom = 4 

f.multifreq = (218, 240, 262, 280, 350) 
f.analyze 0 

RESULTS ["PNA-4 CD 2"] = f.Resuits 
Fig("PNA-4 CD") 

PNA 4 curve could not be fitted 

PNA 4 Tm: 80.2C +/- 1.70 beta: 0.04 cbi2:0.731 

Figure S28 : PNA-4 CD 
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In [24]: f.meth = "CD" 
f.tampon = 1 
f.cell = 4 
f.nom = 6 

f.multifreq = (218, 240, 262, 280, 350) 
f.analyze 0 

RESULTS ["PNA-6 CD 2"] = f.Results 
Fig("PNA-6 CD") 

PNA 6 curve could not be fitted 

PNA 6 Tm: 80.9C +/- 3.60 beta: 0.03 chi2:0.996 

Figure S29 : PNA-6 CD 
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4.3 PNA 5 

Buffer is in cell 1 and PNA-5 in cell 3 

PNA-5 is non-chiral, so CD is irrelevant 
wavelengths : 200 207 215 250 

In [25] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/" 
f.fnom = "PNA_120919_ThermDenat-l-" 
f.extension = ".csv" 
f .IsO 

List of files : 

PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-200nm-cell 1. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-200nm-cell 3. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-207nm-cell 1. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-207nm-cell 3. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-215nm-cell 1. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-215nm-cell 3. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-250nm-cell 1. csv 
PNA_DATA/Dichro_PNA5/Abs_Te/ThermalDenat_120919/PNA_120919_ThermDenat-l-250nm-cell 3. csv 

In [26]: f.meth = "ABS" 
f.tampon = 1 
f.cell = 3 
f.nom = 5 

f.multifreq = (200, 207, 215, 250) 
f.analyze 0 

RESULTS [''PNA-5 Abs"] =f . Results 
Fig ("PNA-5 Absorbance") 
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PNA 5 curve could not be fitted 

PNA 5 Tm: 73.5C +/- 0.54 beta: 0.04 chi2:0.148 

Figure S30 : PNA-5 Absorbance 
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In [27] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_PNA5/Abs_Te/ThermalCooling_120920/" 
f.fnom = "PNA_120920_ThermCool-l-" 
f.extension = ".csv" 

#/. IsO 

f.meth = "ABS" 
f.tampon = 1 
f.cell = 3 
f.nom = 5 

f.multifreq = (200, 207, 215, 250) 
f.analyze 0 

RESULTS ["PNA-5 Abs 2"] =f.Results 
Fig ("PNA-5 Absorbance") 

PNA 5 curve could not be fitted 

PNA 5 Tm: 70.8C +/- 1.77 beta: 0.03 chi2:0.317 

Figure S31 : PNA-5 Absorbance 
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4.4 PNA 7 

PNA 7 is quite special, it presents a very high, much over than the highest experimental temperature 
tested. In consequence, the experimental curves only sample a small portion of the melting curve, and do 
not permit a precise evaluation of the melting temperature. 

It can be seen in the curves, and the results given for the different experiment are not coherent. 

In [28] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130711" 
f.fnom = "PNA8xx_Denat_130711-l-" 
f .IsO 

f.meth = "ABS" 
f.tampon = 1 
f.cell = 5 
f.nom = 7 

f.multifreq = (218, 240, 262, 280) #, 350 ) 

f.analyze 0 

RESULTS ["PNA-7 Abs"]=f .Results 
Fig("PNA-7 Absorbance") 

List of files : 

PNAJ)ATA/Dichro_130711/PNA8xxJDenat_130711-l-218nm-Cell l.txt 
PNA_DATA/Dichro.l30711/PNA8xx_Denat_130711-l-218nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-218nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell l.txt 
PNAJ)ATA/Dichro_130711/PNA8xxJDenat_130711-l-240nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-240nm-Cell 5.txt 
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PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 3.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat.l30711-l-262nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-262nm-Cell 5.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-280nm-Cell 5.txt 
PNAJ)ATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell l.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 2.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 3.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 4.txt 
PNA_DATA/Dichro_130711/PNA8xx_Denat_130711-l-350nm-Cell 5.txt 
PNA 7 curve could not be fitted 

PNA 7 Tm: 90.2C +/- 5.25 beta: 0.16 chi2:0.053 
Figure S32 : PNA-7 Absorbance 



In [29]: f.meth = "CD" 

f.multifreq = (218, 240, 262, 280) 
f.analyze(Tm=80.0,beta=0.05) 
Fig("PNA-7 CD") 

Figure S33 : PNA-7 CD 
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Here, the program is not able to fit the data, even though the initial values (plotted above) seem pretty 
good. 

Probably the value is very high, and cannot be determined. 

In [30] : f = MultiexpO 

f.BASE = "PNA_DATA/Dichro_130711" 
f.fnom = "PNA8xx_Cooling_130711-l-" 

#/. IsO 

f.meth = "ABS" 
f.tampon = 1 
f.cell = 5 
f.nom = 7 

f.multifreq = (218, 240, 262, 280) 

f.analyze 0 

RESULTS [''PNA-7 Abs 2"] =f .Results 
Fig("PNA-7 Absorbance”) 

PNA 7 curve could not be fitted 

PNA 7 Tm: 72.1C +/- 3.24 beta: 0.15 chi2:0.064 
Figure S34 : PNA-7 Absorbance 
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In [31]: f.meth = "CD" 

f.multifreq = (218, 240, 262, 280, 350) 

f.analyze 0 

RESULTS ["PNA-7 CD 2"] =f.Results 
Fig("PNA-7 CD") 

PNA 7 curve could not be fitted 

PNA 7 Tm: 165.5C +/- 149.85 beta: 0.02 chi2:0.501 
Figure S35 : PNA-7 CD 
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5 conclusion 

Results have been accumulated and are given below: 

In [32] : n^print RESULTS 

for k,v in RESULTS . items(): 

print "yoS Tm : ydf +/- yo.2f"yo(k, v[^Tm^], v [ ^Tm_errorbar ^ ]) 


PNA-2 

Abs Tm : 

66.3 +/- 1.30 

PNA-3 

Abs Tm : 

71.1 +/- 1.77 

PNA-3 

CD Tm : 

68.4 +/- 0.98 

PNA-3 

CD 2 Tm 

: 68.1 +/- 1.82 

PNA-4 

Abs Tm : 

86.4 +/- 4.18 

PNA-6 

Abs Tm : 

88.8 +/- 5.51 

PNA-4 

CD Tm : 

81.7 +/- 2.43 

PNA-6 

CD Tm : 

83.0 +/- 4.50 

PNA-4 

CD 2 Tm 

: 80.2 +/- 1.70 

PNA-6 

CD 2 Tm 

: 80.9 +/- 3.60 

PNA-5 

Abs Tm : 

73.5 +/- 0.54 

PNA-5 

Abs 2 Tm 

: 70.8 +/- 1.77 

PNA-7 

Abs Tm : 

90.2 +/- 5.25 

PNA-7 

Abs 2 Tm 

: 72.1 +/- 3.24 

PNA-7 

CD 2 Tm 

: 165.5 +/- 149.85 
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